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This paper presents a simplified set of equations 
for calculating and optimizing regenerator geometries. 

A number of competing parameters oust be accounted for 
in the design of a regenerator. To obtain high values 
of effectiveness, there is need for a large surface 
area for heat transfer and high heat capacity. The void 
volume should be small to maintain the pressure ratio in 
the entire system. Moreover, the pressure drop should 
be small compared to the absolute pressure. With the 
assumptions made here, the calculations can be done with 
a hand calculator. 

The equations show that the optimum regenerator length, 
hydraulic radius, and porosity are independent of mass flow 
rate and the gas cross-sectional area is proportional to 
the mass flow rate. It is shown that using gas gaps 
between parallel plates produces a significantly better 
regenerator than is possible with packed spheres or screens, 
particularly at temperatures below approximately 50 K. 

The improvement is due to enhanced heat transfer and the 
flexibility in using a lower porosity in the gap configur- 
ation. ' 
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1. Introduction 


The regenerator is often one of the major loss sources in regenerative-cycle refrigerators. 
It contributes loss terms due to its limited heat transfer units, limited matrix specific heat, 
pressure' drop, dead volume, and axial thermal conduction. An optimum design of these re- 
generators would significantly improve the performance of the overall refrigerator. 

The calculation of regenerator performance with various assumptions has been discussed exten 
sively in the literature. First order calculations assume the void volume in the regenerator is 
zero. Such calculations were first done by Lamberts on [1] and Hausen [2] with graphs and tables 
presented in the books by Kays and London [31 and Schmidt and Willmott [4]. Second order calcu- 
lations take into account the regenerator void volume but neglect any pressure oscillation in the 
void volume gas. These calculations have been done by Heggs and Carpenter [5], by Daney and 
Radebaugh [6], and by others cited in these two works. A third order calculation considers void 
volume in the regenerator and a pressure oscillation as seen by regenerators in regenerative- 
cycle refrigerators. Gary, et al . [7] discuss such a third order calculation where very few as- 
sumptions are made. 


*Work partially funded by the David Taylor Naval Ship R&D Center and the Office of Naval 
Research. 
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All of the. cited calculations require the regenerator geometry as the input parameter and 
compute the regenerator effectiveness. The problem facing the designer is how to determine the 
optimum geometry without resorting to trial and error through such performance calculations. 

Such an approach can be very costly and time consuming, especially for third order calculations 
where considerable computer time js required for each case. In this paper we discuss a simple 
technique for determining the optimum geometry associated with a desired value of regenerator 
effectiveness. The technique reverses the "normal" first order calculations and computes the 
various geometrical parameters that give the minimum regenerator void volume.. If the void volume 
is not negligible a second iteration should be done with reversed second or third order cal cula- 



2. Basis for optimization 


In any optimization scheme it is important to determine the proper parameter to be optir. 
mized. Most often the input power is minimized for a refrigerator with a specified net refriger- 
ation power. (In some applications, such as cooling SQUID'S, the net refrigeration power is zero 
and the gross refrigeration is used entirely for loss terms. Nevertheless, the input power is 
still the quantity to be minimized in SQUIO cryocoolers. ) It would appear that a general system 
requirement would preclude a simple optimization of a single component like a regenerator, espe- 
cially when the regenerator may consist of several stages. Any exact analysis must consider the 
entire system as a coupled system and various loss terms will necessarily be dependent on each 
other. However, Smith and coworkers [8-11] have found that a decoupled approach to the loss terms 
has resulted in relatively good approximations to the overall system performance. In the decoup- 
led approximation the gross refrigeration power, 0 r » is absorbed by several independent terms: 


0 = Q + 0 . + 6 + 6 + 6+6, 

x r m, 'rad y c x reg N h 


( 1 ) 


where Q net is the net refrigeration power, Q pad is the radiation loss, Q c is the conduction loss, 
Q re g is the loss due to regenerator ineffectiveness, Q s is the shuttle heat loss due to the oscil- 
lating displacer, and Q h is the loss caused by an excess enthalpy flow through the regenerator 
during the hot blow for certain conditions with a non-ideal gas. Another possible term would be 
one due to the pressure drop, AP, in the system. That term would be used whenever Q p is 
calculated from the pressure amplitude at the compressor. If Q p is calculated from the pressure 
amplitude at the expansion space, then no such term would be used, but the pressure drop would 
indicate the required pressure amplitude at the compressor. 

In the decoupled approximation of eq.(l) the entire system is optimized by minimizing Q p for 

a fixed value of Q net - For an ideal gas Q^ = 0. For the second and third stages 6 raci is often 

negligible. The other terms, Q , <5 , and Q , are then usually of comparable vali ;s when a min- 

- , ** ' s 

imum Q p is found. Likewise, the percentage loss due to a pressure drop, AP/P, will be comparable 
to the percentage loss of these other three terms. The simplest approximation for 6 r assumes an 
isothermal expansion. If an adiabatic expansion is considered, Q p is reduced by approximately a 
factor of two [12]. 
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In dimensionless form eq.(l) is given as 


*neA 


<W*r + Mr + fi re A 


0 /Q 


Mr 


= 1 


( 2 ) 


For any particular refrigerator an experienced designer can then make reasonable estimates of 
such terms as Q c /Q r and Q re g/Q r and AP/P. For example, with an ideal gas in a cryocooler where 
^net^r to be sma ^» t * le aforementioned terms may be approximately 0.1-0. 2. Good estimates of 
these terms are the first and most important step in the optimization of a regenerator. The per- 
formance of the overall system is also affected by the void volume of the regenerator, V rg . A 
normalized void volume would be V r g/V g , where V g is the expansion space volume. The four terms 
Q r eg/Q r > Q c ^r* ap/p > antl V rg^ V e are interre1ated ‘> a decrease of one must cause an increase of 
another. 

The basis for the optimization procedure described herein is to use O re g/Q r » 0 c /$ r . and &p / p 
as known input parameters and then find the geometry which gives a minimum V /V fi . At an early, 
intermediate point, the regenerator ineffectiveness assoicated with 0 reg /§ r 1S calculated. If 
V rg /V e « 1, then the input parameters could be decreased to the point where v r g/ v e becomes 
significant. If V rg /V fi >> 1. then the three input parameters must be increased. ' If the values 
are as large as they can be.^the process probably is not feasible with the conditions chosen, 
unless an abnormally large compressor is used. The optimization done here is not fully rigorous 
since estimates are used for the input parameters and the effects of it P/P and v rg / v e on the 
system are not fully quantified. 

The geometry of a regenerator can be divided into two components: (a) the configuration of 

the packing, such as packed spheres, plates, tubes, etc. and (b) the macroscopic geometry, such 
as cross-sectional area, porosity, length, and hydraulic radius. The optimization procedure 
discussed here treats both components. The best configuration is the one which has the highest 
heat transfer rate for a fixed pressure drop. However, any configuration may be chosen and the 
optimization procedure then gives the best. geometrical parameters. Pressure drop and heat 
transfer correlations often incorporate the friction factor, f, and the Stanton number 


H st = h/(m/A g )c p ; 


(3) 


where h is the heat transfer coefficient, m is the mass flow rate, A g is the gas cross-sectional 
area, and c p is the specific heat of the gas at constant pressure. 

2/3 

Figure 1 shows a plot of N N and f as a function of the Reynolds number, N , for one 

si. pr o /3 

particular configuration, an infinitely long and wide gap. The two curves of N st N p ^. are for the 

case of constant heat flow and for constant temperature. Kays and London [3] provide the 

correlations for gaps and several other configurations. The Prandtl number N pr is included to 

make the correlations valid for a wide range of fluids. For helium gas N = 0.666 under ideal 

2/3 pr 

conditions. Note that both and f have nearly the same dependence on N r - For 

regenerators the curve of interest is that of constant heat flow. The practical difficulty in 
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achieving a uniform hydraulic radius, however, suggests the lower, constant temperature curve Bay 
he more realistic. ' 


We define the ratio: 


• ■ V# 3 " 


(4) 


A high value of a is desired since it gives a high heat transfer for a given pressure drop. (A 
more quantitative argument is given later.) Figure 2 shows a as a function of the Reynolds num- 
ber for several packing configurations and for a constant temperature heat transfer. The data 
are from Kays and London [3]. The two important points to note from these curves are fa) they 
are nearly independent of Reynolds number, and (b) a gap configuration has the highest or - sig- 
nificantly higher than that of a packed bed of spheres. Because a is nearly independent of 
Reynolds number, it becomes a useful parameter to use in calculations since there is no need to 
specify the flow rate. 



1(f 


Figure 1. Heat transfer and friction factor Figure 2. Ratio of Heat transfer and friction 
curves as a function of Reynolds number for factor curves as a function of Reynolds number 
gas flow in a gap. for several configurations. 


3. Optimization procedure 


The equations which are developed in the following sections describe the procedure to calcu- 
late the actual dimensions of an optimum regenerator geometry. The development begins by relating 
Q c /Q r to a regenerator ineffectiveness. By using published performance curves for regenerators, 
a set of values for the number of heat transfer units per half cycle, N^ u , and the ratio of matrix 
heat capacity to the heat capacity of the fluid that passes through the regenerator, C r /C f , can 
be found that yields the correct ineffectiveness. Because an infinite set of values will yield 
the same ineffectiveness, the goal is to find the set of and C r /C^ values that gives the min- 
imum V_ /V . Expressions for the geometrical parameters of cross-sectional area, length, 
r g e . 

porosity, and hydraulic radius are developed in terms of the optimized N t and C r /C f and the 

Input parameters Q c /Q r and AP/P. The configurational parameter a is used in the equations to 

make them valid for all N and to show how a maximum or produces a minimum V /V . 

r r rg e 
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N tu 


Figure 3. Calculated regenerator ineffec- Figure 4. Calculated regenerator ineffec- 
tiveness for zero void volume. tiveness for a finite void volume but with no 

pressure wave. 

3. 1 Performance curves 

The effectiveness, e, of a regenerator is defined as the ratio of heat actually transferred 
to the maximum possible heat transfer. The first optimization step utilizes the performance 
curves as shown in figures 3 and 4. Figure 3 shows the ineffectiveness, 1 - e, as a function of 
N tu and C r /C^. The plot assumes zero void volume, zero thermal conductivity in the axial 
direction, infinite thermal conductivity in the radial direction, gas and matrix thermal 
properties independent of temperature, constant inlet temperatures, constant flow rate, and 
constant heat transfer coefficient. The curves in figure 3 are plotted from data in Kays and 
London [3] modified by. the appropriate change for our definition of N t(j . Figure 4 is for the 
case of a finite void volume without a pressure wave [5,6]. All other assumptions are the sasse 
as for figure 3. Other curves could be generated using fewer assumptions. Since the thermal 
conductivity in the axial direction is taken into account separately via a decoupled loss terra, 
it should always be neglected in the performance curves. 

For a predetermined value of O re g/Q r a value of 1 - e can be calculated for a regenerator 
operating between some upper temperature T^ and a lower temperature T^. By definition 

1 ‘ 6 = W“ (H o ' H L } * (5) 

where H is the enthalpy per unit mass. Since refrigeration in a regenerative-cycle cryocooler 
occurs only during one half-cycle, the refrigeration rate is given as 

Q r = ™l r /2, (6) 
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where q p is the heat absorbed per unit mass during the expansion process. Combining eqs. (5) and 
(6) gives 

1 ' * = (Qreg/VV^U ' V\ (7) 

As shown in figure 3, a given ineffectiveness is satisfied by a wide range of and C r /C^ 
values. The optimum combination will become evident later, but we proceed with the optimization 
procedure assuming that the best N^ u and C r /C^ have already been determined. 


3.2 Gas cross-sectional area, A 

9 

A high in a heat exchanger is achieved by using a high surface area, but the high 
surface area will also lead to a large pressure drop unless the cross-sectional area, is 
sufficiently large. The pressure drop in the core of the regenerator is given by 


AP = f(m/A g ) 2 L/2pr h , 


(8) 

where f is the friction factor, L is the regenerator length, p is the gas density, and r^ is the 
hydraulic radius. Equation (8) does not consider the entrance and exit pressure changes associ- 
ated with acceleration and deceleration, but these terms are generally negligible in a well de- 
signed regenerator. By using the term a from eq. (4) the pressure drop in eq. (tj can be related 
to the heat transfer by 

* P = N st N pr 3(,;,/ y 2 L/2of P r h- . 


(9) 

By using the definition for in eq. (3) and the following definitions 



N. = hA/rnc 

VsJ P 


(10) 

r h =L yV 


(11) 


where A is the total heat transfer area, we can express N^ u in terms of N gt as 


N tu = ;y, st L/ V 


( 12 ) 
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Rearranging eq. (12) and substituting into eq. (9) for H Jt allows us to write for the gas 
cross-sectional area 


V“ = [N tu N^ 3 /2 a pAP) !i . 


( 13 ) 


Equation (13) is one, of the most important equations in the design of a heat exchanger, whether 
it is regenerative or recuperative. It gives the gas cross-sectional area needed to achieve a 
given for a fixed value of AP across the exchanger. Equation (13) also shows that should 
be linearly proportional to ro and that a large o means a small A g . The term a is nearly constant 
for a given heat exchanger packing and is taken from figure 2. An average value of p must be 
used when large temperature differences occur from one end to the other. In dealing with a mass 
flow rate that is r.ot constant, e.g. sinusoidal flow, the average m should be used to determine 
Ag. Likewise the AP then corresponds to the average m and not the peak. 

3.3 Regenerator length, L 


Equation (13) shows that V ro could be made as small as desired by using a short regenerator. 
However conduction in the regenerator would then become large. We now calculate the length 
assoicated with the given value of Q c /Q r - 8y definition 



(14) 


where A^ is the regenerator solid or matrix cross-sectional area and k is its thermal conductiv- 
ity. The term is related to through the porosity n g by 


n = A /(A 
9 9 ■ ra 


A ). 
9 


- (15) 


Equations (15) and, (6) are substituted into cq. (14) to yield 

i 

2(A /m)(l - n ) JkdT 

L = 3 — -S . ( 1 6) 

n q CO /Q ) 
g M r' x c H r' 

For some packing configurations is fixed, such as packed spheres or stacked screens. In those 
cases, L can then be calculated from the known input parameter Q c /Q r - For other configurations, 
such as gaps, the porosity is entirely flexible and L cannot be calculated until the porosity is 
first calculated. Hote that eq. (16) also applies to either recuperative or regenerative heat 
exchangers, since the term has not been used. 
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3.4 Porosity, n 

0 

In configurations where the porosity is flexible, additional matrix heat capacity can be 
obtained by adding more matrix material and still keeping A g constant. However, at some point 
conduction effects become important. The regenerator matrix heat capacity is given by 


C 


r 


Vm c m* 


( 17 ) 


where p is the matrix density, c 
m m 


is the matrix specific heat, and V is the matrix volume 


V = 
m 


At . 
m 


(18) 


Equation (15) is used to express 


in terms of A 

9 


by 


V m = A L(1 - n )/n . 
mg g 9 


(19) 


Substitution of eq. (19) into eq. (17) gives 


C r = 


<VmV (1 


n )/n . 
9 9 


( 20 ) 


The heat capacity of the fluid that passes through the regenerator is 


= mCpt/2 = mc p /2v, (21) 

where x is the period of one cycle and u is the frequency of operation for the regenerative cryo- 
cooler. The use of eq. (21) now restricts subsequent results to regenerative heat exchangers. 
When eq. (20) is divided by eq. (21) and rearranged the result is 

n g /(i - n g ) = 2u(p ra c m )(A g /m)L/c p (C r /C f ). (22) 


For simplicity in subsequent calculations we define X as 
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X - n g /(l-n g ) 


n g = X/(l ♦ X) 


The length in eq. (16) is substituted into eq. (22) to yield 


X = 2(A Ja) 


v( Pm c m^ kdT 

c p q r ( c r / C f )(Q c /Q r ) 


In calculating n g the ratio C r /C^ is assumed to be known froiu the performance curves in figure 3, 
although any combination of ll tij and C r /Cf can be used. The optimum combination is yet to be 
determined. Equations (24) and (25) show that the porosity is independent of the refrigerator 
capacity since (A g /m) and the other terms are independent of mass flow rate. Likewise the length 
in eq. (16) is independent of mass flow rate or refrigerator capacity. Also, the porosity 
equation applies to the case of variable porosity. The proper technique for dealing with a fixed 
porosity is discussed in section 5. 


3.5 Regenerator gas volume, V pg 

With A and L already determined, V is easily calculated. However the ratio V /V needs to be 

g rg J rg e 

calculated to determine what effect V pg will have on the system. The volume of gas in the 
regenerator void space can be given by 


V = XV . 
rg m 


The matrix volume can be expressed as 


^r' /p m C m’ 


which means eq. (26) becomes 


V = XC /p c . 
rg r m 
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The expansion space volume is given by 


V e = CfVp* 


(29) 


where p e is the density of gas in the expansion space. The tern p should be evaluated at the 
temperature and pressure of the expansion space at the time of flow reversal. For a first approx- 
imation the lower temperature T^ and the average pressure are used. From eqs. (28) and (29) the 
volume ratio becomes 


V 

rg 


/V 


e 




(30) 


This simple expression for V^/V^ is useful when X has already been calculated. However, since 
the optimum N^ u and C r /C^ have not been selected, we now develop an expression, for v r g/ v e in 
terms of N t and C r /C f . We substitute eq. (25) for X into eq. (30), and at the same time eq. 
(13) is used for A^/m in the expression for X. The resultant volume ratio becomes 


V /V = 
rg e e 


2N tu< C r /C f>Vr %J kdT 

apAPq r (Q c /$ r )(P m c ffl ) 


(31) 


It now becomes clear how V r g/V g depends on the various input and material parameters and the 
packing configuration. The volume ratio is minimized when the configurations! parameter or is 
largest. Figure 2 shows that gaps offer the largest a for the conf igurations compared. The 
authors are unaware of other configurations with a higher or. Of course, the equations derived 
here are valid if one chooses to use a configuration with a smaller or. In regard to matrix 

material parameters, the ratio [JkdT/(p m c )] should be minimized. The input parameters 6P/P "and 

0 C /Q r should be made large to reduce V r gA' e but they cannot be made larger than about 0.1-0. 2. 

A minimum value for V /y e is obtained when the product N tlf (C r /C^) is a minimum. Figure 5 

shows curves of N^ u (C p /C^) vs. (C^/C^) for various values of ineffectiveness. We see from figure 
5 that C r /Cf = 1.8 gives a minimum in N t (C /Cf). for all values of 1 - c shown. The curves of 
figure 5 are subject to some uncertainty since they were derived from the graph in figure 3 rather 
than using numerical calculations. Fortunately the minimum in N tu ( c r / C f) is fairly broad so some 
uncertainty can be tolerated. The calculation of (A^/ro), L, n^, and V p g/V e is now done using the 
optimum and C p /C^ combination from figure 5. 

For an ideal gas and isothermal expansion we have 

P = P/RT i 32) 
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Figure 5. The N tu (C r /C r ) product for several values of regenerator ineffectiveness 
from the case of zero void volume. 


P D = P/RT„ 
e e 


(33) 


q r = RT e £n(P h /P £ ), (34) 

where n is tne gas constant, P is the average pressure, P^ is the high pressure, P £ is the low 
pressure, and is the' expansion space temperature. (Usually T e = T^). Then eg. (31) becomes 


V /V = 
rg e 



2N tu (C r /C f^ N pr S T ^ kdT 

V (p mSs )( VV (aP/P)1n(p h /, V 


(35) 


for an ideal gas and isothermal expansion. 

Equation (35) shows that the volume ratio is independent of pressure and 
will increase as the temperature T g is lowered. For the optimum (C^/C^) = 1.8, the 
ineffectiveness is approximately 

1 * e = 2. 1N~®' 9 , (36) 
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f°r_N tu > 50. The tern N ty is then proportional to (1 - e)’*‘* and using eq. (7) for (1 - e) 

shows N. is proportional to (6 a „/6) The denominator inside the brackets of eqs. (31) and 

tu rsy r • • 

(35) subsequently contains the three input parameters in the forta of (0.„„/$„) i ‘ X (Q*/0-)(dP/P). 

icy r c r 

Such an expression would suggest that ($ re g/$ r ) should be larger than the other terns for mini- 
mizing v rg /v e . 

There may be instances where the total regenerator volume is to be minimized. This possibil- 
ity may occur when there are external packaging constraints or in cases where no pressure wave 
exists. For instance, regenerators used in magnetic refrigerators operate at constant pressure. 
Here, the void volume is often of little concern. Of more, importance is the total regenerator 
volume, since it may need to fit inside a magnetic field region. Tho total regenerator volume, 

V », can be expressed in reduced units by 


V ./V = V /V n . 
rt e rg e g 


(37) 


The use of eqs. (24), (25), and (30) leads to 


P e (C p /C f )c 

V ./V = ■ / ■■ -■ v P - + p 

rt e < p m c J e 


Wf^ 3 ^ 


«pAPq r (Q c /O r )(P m C n ) 


(38) 


The first term on the right hand side of eq. (38) is the reduced matrix volume and the second 
term is V rg /V e of (31). It is now apparent that the minimum in V rt /V e occurs at a smaller 
value of (C r /C f ) than the 1.8 used for V rg /V e , and the f« tu will be higher than for that case. 


3,6 Hydraulic radius 

Here we derive the hydraulic radius and other characteristic dimensions of the regenerator 

packing that must be used to give the desired in a regenerator with a length and 

cross-sectional area as previously calculated. When both sides of eq. (12> are multiplied by 
2/3 

N and rearranged, the result is 


" s t»pr 3 ’ VprV 


(39) 


The Reynolds number is defined as 


N r = 4r h /p(A g /m), 


(40) 
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wfjere p Is the gas viscosity. 

For every geometry of interest there is some function g where 


"stV = S (N r>' {41 > 

2/3 

Equations (39-41) represent three equations in the three unknowns, r^, fi^, and M st Np r that ®ust 

be solved numerically in the generat case. When g is in graphical form, the solution is easily 

2/3 

found by trial and error. For a chosen value of r. , eq. (39) is used to find N *N_ _ . The same 

n 2/3 P" 

value of r. gives H from eq. (40) which in turn gives a value of M from the graph of eq. 

n r 0/3 St pr 

(41). If the two values of K do not agree, another value of r. is chosen until the two 

2/3 st pr h 

N st N pr Wfilues a 9 ree - 

For laminar flow , eq. (41) is 


N st N P ? /3 = a/N r* < 42 > 

where a is a constant depending on geometry. For infinitely wide and long gaps with constant 
temperature heat transfer, a = 8.5. For infinitely long tubes, a = 4.2. Equations (39), (40), 
and (41) are solved then by 


r h = 


alp(A /m) 

VT~ 

4N* fi J 
tu pr 


(43) 


In the high-performance heat exchangers discussed here, the flow is usually laminar. The value 
of or in figure 2 can be refined with if necessary. Since a is only a weak function of H , 
further iterations may not be necessary. 


The characteristic dimension of a heat exchanger may be expressed in terms of the hydraulic 
radius. They are given for various geometries as follows: 


tg = 2r h gap thickness (44) 

d = 4r h tuba diameter (45) 
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d = 6i* h /X sphere diameter (46) 

d = 4r f) /X wire diameter. (47) 

for a gap configuration the thickness of the plates that comprise the matrix is 

t = t/X. (48) 

m g 

For all of the preceding calculations to be valid, the thermal penetration depth X must be 
large enough to penetrate throughout the matrix. For a serai-infinite plate the theraal penetra- 
tion depth for sinusoidal heat flow is 

A = tk/<P B c a )n»]\ (49) 

For a set of stacked plates with gas on both sides we require 

t < 2k. (50) 

ID 


3.7 Equation summary 


As an aid to computations for regenerators, the significant steps and equations are summarized in 
the order necessary for a regenerator with flexible porosity. 

1. Determine largest input parameters 0 re g/§ r » § c /§ r . the overall system can tole- 
rate. 


2. Evaluate system and material parameters v, Tg, T^, 'a, >\ q r> p, p e , c p , N pr , p, (P ra c m ), 
and JkdT. 

3. Ineffectiveness: 1 - £ = * H,_). 

4. From figure 3 or 5 determine: if. at C /C~ =1.8. 

tu r T 


5. Determine packing parameters a from figure 2 and laminar coefficient a in eq. (42). 

2/3 k 

6. Gas cross-sectional area: A /m = [M. N /2apAP] . 

g tu pr 


7. Porosity: n g = X/(l + X) 

X = 2(A /ra) 


' u( Pfil c n );kd T 

c pV C r /C f )( W 
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(in original parameters) X = 2 




2opAPc p q r (C r /C f )(Q c /Q r ) j 
8. Regenerator gas volume: V rg /V g = X(C r /C f )c p P e /(P m c o ). 


(in original parameters) 


V /V = p 
rg e 


Z \^ /C yjrW MX 


1 


apAPq r ($ c /q r )(P m /c n ) 


9; Regenerator length: L = 2(A g /m)J'kdT/Xq [ ,(Q c /() r ). 

I s * 


(in original parameters) 


L = 


c p (C r /C f )/icdT 


u(p c )q (Q /Q ) 


faLp(A /m) j 5 

10. Hydraulic radius: r h = ] — 


L4N tuV 


laminar flow. 


(in original parameters) 


r h = (ap)' 


c p (c P /C f W kdT 


1 h 




11. Reynolds number: N p = 4r h /p(A g /m) 

Recalculate r h if flow not laminar. 

12. Reevaluate a. If necessary, repeat from step 6. 

13. For gaps: t = 2r. , t = t /X 

31 g h m g 

s 

14. Thermal penetration depth in plates: K = [k/(p m c m )nv) ' 

sinusoidal flow. 


15. Check to see if: t < 2A. 

m 

If not, increase u to decrease t . 

m 

16. Error check: (a) AP = fi./Zpr^(A g /m) 2 , 

with f from figure 1 or similar curve. 

' Calculated value should equal input value. 

(b)V rg /V e = 2(A g /m)Lup e , 
should be same as from step 8. 


( 51 ) 


(52) 


(53) 


4. Example solutions 

The solution to the design equations can best be illustrated through some examples. In the 

design of a cryocooler for a SQUID, nonmagnetic and nonraetallic parts need to be used in the cold 

parts to prevent magnetic noise. The regenerator material chosen for this illustration is 6-10 

fiberglass epoxy. Even though its volumetric heat capacity is small compared with lead and other 

regenerator materials, its very lew thermal conductivity gives it a desirably small value of 

k/p^c^, the important material parameter which appears in eq. (31). In fact G-10 may have a 

lower k/p c than Pb-5% Sb although k for the latter has not been measured, 
m m 
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The parameters '•hosen for the example solution are given in table 1. 


Table 1. Input parameters for example solutions 


1 - e = 0.01 

Q c /Q r = 0.1 
AP/P = 0.05 
P = 0.5 MPa 
P h /P £ = 2.0 


(Q /Q = 0.22) 
vv reg y r ' 



or = 0. 35 (gap) 
regenerator material: G-10 

fluid: ideal helium gas 

isothermal expansion 


The temperature T^ is allowed to vary from 4 K to 100 K to evaluate the temperature dependence of 

the optimum geometry. The upper temperature T^ of the regenerator is assumed to be 4T^, which is 

typical for cryocoolers. Such a dependence of T^ on T L means that eq. (7) may be written as 

1 - e = <W*rV 6c p T L (54) 

for the case of an ideal gas. Also, for an ideal gas q f from eq. (34) can be substituted into 
eq. (54), and for the normal case of T g = T^, we obtain 

1 - c = «reg^r )R£n(P h /P P /6c p- (55) 

For the input parameters in table 1, (Q rg g/Q r ) = 0.22 from eq. (55). 

The material parameters needed in the calculations are (pc) and flcdT. Figure 6 shows 

m nr 

(p R c ra ) for G-10 fiberglass-epoxy, Pb [13], and GdRh [14]. The alloy Pb + 5% Sb is commonly used 
as a regenerator material and GdRh is a potential regenerator materia) because of its high speci- 
fic heat (P m c m ) at low temperatures. The thermal conductivity of G-10 is the average of the paral 
lei and perpendicular heat- flow directions from the work of Kasen, et al. [15]. Table 2 gives 
representative values of (p m c m ), k, and |kdT used in our calculations. 

For the ideal helium gas we use c p = 5.19 J/g-K, R = 2.077 J/g-K, N pr = 0.665. The viscosity 
values from McCarty [16] are nearly independent of pressure over the range we have considered and 
are approximated by p = 5.0 x 10 ® + 4.6 x 10 ®T^'^ g/cm.s over the range of 10-300 K. 

Figure 7 shows (A^/m) and the 1-ength of the regenerator as a function of T^ for three differ- 
ent frequencies. The term (A^/m) is independent of frequency. The porosity, n^, and the volume 

ratio, V /V , are shown in figure 8. 
rg e’ J 

Note that the optimum porosity is on the order of 0.01 for T^ = 5K. Such a low porosity 
allows a large volume of regenerator material to be used to obtain a reasonably high total heat 
capacity. The volume ratio V /V is insignificant at higher temperatures and remains smaller 

*9 
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TEMPERATURE. K 


Figure 6. The volumetric heat capacity of candidate regenerator materials. 


Table 2. Thermal properties of G-10 fiberglass epoxy 


T 

(K) 

P c 
K m m 

(J/cm 3 -K) 

k 

(W/cm-K) 

f 0 kdT 

(Vl/cm) 

5 

0.0058 

0.815xl0' 3 

2. ODxlO* 3 

10 

0.0255 

1.24 xlO* 3 

7.23X10" 3 

20 

0.078 

1.73 xlO* 3 

22.3 xlO* 3 

50 

0.245 

2.65 xlO" 3 

89.1 xlO* 3 

100 

0.520 

3.76 xlO" 3 

251. xlO* 3 

200 

1.10 

5.53 xlO* 3 

718. xlO* 3 

300 

1.64 

7.17 xlO* 3 

1355. xlO* 3 


than 1.0 even at T^ = 5K. For the case of u = 1 Hz the volume ratio is small enough that the 
assumption of zero void volume in the performance curve of figure 3 may hold. For higher fre- 
quencies, a more realistic performance curve probably should be used for the lowest temperatures 
The small values of V r g/V e a ^ s0 imply that the entire process for refrigeration at 4-5 K is feas 
ible with an ideal gas. (Calculations for a real gas are easily done but are not shown here.) 

Figure 9 shows the gap and matrix thicknesses as a function of T^. Shown for comparison 
with the matrix thickness, t , is the thermal penetration depth, 2\. For all cases except at 1 
Hz and the lowest temperature << 2 K. When t s 2X, there will be a slight degredation in re- 
generator performance because the term C will be reduced. The gap thickness has a slight deperi 

' -3 

dence on temperature and is in the range of 10 cm. An extremely small t could present some 
practical problems in the construction of the regenerator. The Reynolds number also is very 
small - in the range of 10-200 - with higher values occurring at lower temperatures and lower 
frequencies. 
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Ag/m, cm 2 s/g 


21 


100 
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5. Comparison of gap and packed sphere configurations 

Packed spheres are coestonly used for regenerator packings since such regenerators are sl^sle 
to make. However a is significantly lower for packed spheres than for gaps. Also, the porosity 
of packed spheres is fixed at about 0.38, which is isuch higher than optimum for temperatures be- 
low about 50 K according tc figure 8. In this section we show how auph inferior a packed sphere 
bed is in comparison with the gap configuration. Of course, the packed sphere bed still has the 
advantage of simplicity. For a fixed n , eq. (51) is solved for the ratio N tu /(C r /Cf). The 
optimum values for N t and C r /C^ are found by simultaneous solution with the curves of 

figure 5. This approach maintains (Q c /Q r ) at the fixed input of value. Included in the 
assumptions is that /fcdT for spheres is the same as for the bulk material. After and C f /C^ 
are determined, the parameters A^/ra, L, V /V e , and r^ are calculated from the appropriate 
equations. 


Tu.* 

10 20 SO 100 200 



Figure 10. The reduced regenerator gas volume for gaps and for packed spheres. See test for 
details of packed-sphere curves. 
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Figure 10 shows V r ^V g * or both packed spheres of G-10 and gaps between plates of G-10. At 
higher temperatures the difference between the two sets of curves is mainly due to the higher a 
for gaps. At lower temperatures the such smaller V r g/V g for gaps is due to both the higher o and 
the lower porosity. Also shown in figure 10 for packed spheres is the curve labeled = 378, 
which is the value used for the gaps. At higher temperatures this curve gives a lower V^/V^ 
because § c /0 r is greater than 0.1. At lower temperatures the reverse is true, and § C /Q r is less 
than 0.1. Wherever the packed-sphere curves are below the = 378 curve, the is greater 

than 378. In fact, for temperatures less than 15 K the fl. is rapidly approaching nusbers as 

5 ^ 

high as 10 . a practical impossibility. Even an N^ u of 378 nay bs difficult to achieve in prac- 
tice without extreme care in making all flow channels the same size. If a saaller JTkdT is used 
in the calculations of the packed powder, an even higher N t(j is needed for the optimum condition. 
However, the reduction in V /V e * s i ns ’9 n ^Fi cant at temperatures below about 15 K. The f« tu = 

378 curve should be considered a lower practical limit for the V r g/V e of the packed spheres. 


6. Conclusions 

A simple scheaie for optimizing regenerators in cryocoolers has been developed which uses the 
system loss terms, Q re g/Q r . Q c /Q p , and AP/P as input parameters. Any set of regenerator perfor- 
mance curves can be used, although the examples discussed here use the simplest one of zero void 
volume. A set of equations ware developed to detereine the regenerator geometry that yields the 
oinieua void votuaa. In comparison -.^ith other packing configurations It is shown that gaps be- 
tween parallel plates is the best, producing regenerator dead volumes one to two orders of magni- 
tude smaller for temperatures in the range of 5-10 K. It is shown that G-10 fiberglass-epoxy is 
a favorable saterial for a multiple gap regenerator because of its relatively low k/p Q c r) ratio. 


The authors are grateful to David Daney arid James Zimmerman, both of KBS, for some valuable dis- 
cussions. Debra Schlender is acknowledged for careful preparation ef the manuscript. 
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8. Nomenclature 


English letter symbols 


SI units 


A 



a 




k 

L 

ra 


regenerator total heat transfer area 

gas- cross-sectional area 

matrix cross-sectional area 

geometrical constant for laminar flow 

heat capacity of fluid passed through regenerator 

heat capacity of regenerator matrix 

heat capacity of gas in regenerator void volume 

specific heat of matrix 

specific heat at constant pressure of gas 

tube, sphere, or wire diameter 

friction factor 

specific enthalpies at T^ and T^ 
heat transfer coefficient 
thermal conductivity of matrix 
length of regenerator 

mass flow rate of gas through regenerator 


J/K 

J/K 

J/K 

J/kg-K 

J/kg-K 

(!) 

J/kg 

W/m 2 -K 

W/m-K 

m 

kg/s 


197 



V 

Prandtl number, a fluid properties modulus 


"r 

Reynolds number, a flow modulus 


N st 

Stanton number, a heat transfer modulus 


N tu 

Number of heat transfer units per half cycle 


"o 

porosity of regenerator = A^/(Ag ♦ A^) 


p 

average pressure 

Pa 

p p 
h* v l 

high and low pressures 

Pa 

AP 

pressure drop in regenerator 

Pa 

% 

conduction heat loss 

W 


enthalpy deficit loss 

W 

<U 

net refrigeration power 

W 

*r 

gross refrigeration power 

W 

^rad 

radiation heat loss 

W 

A reg 

regenerator loss 

w 

$s 

shuttle heat loss 

w 

q r 

heat absorbed per unit mass during expansion 

o/kg 

R 

gas constant 

J/kg-K 

r h 

hydraulic radius (LA^/A) 

m 

T e 

expansion space temperature 

K 

V t l 

upper and lower temperatures of regenerator 

K 

S 

gap thickness 

ra 

* 0 

thickness of plates in gap regenerator 

is 

V e 

maximum expansion space volume 

ra 3 

V m 

matrix volume 

m 3 

% 

regenerator void volume 

B 3 

3 

V rt 

regenerator total volume 

nr 

X 

porosity parameter [n^/(l-n^)] 


Greek 

letter symbols 


a 

2/3 

ratio of heat transfer and friction terms /f) 


€ 

regenerator effectiveness 


1-e 

regenerator ineffectiveness 


A 

thermal penetration depth in matrix 

R) 

M 

gas viscosity 

kg/m-s 

V 

cycle frequency 

Hz 

P 

gas density in regenerator . 

kg/ra 3 

p e 

gas density in expansion space 

kg/s 3 

p m 

matrix density 

kg/s 3 

t 

period of one cycle 

s 
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